Review of measurements of Kerr nonlinearities in lithium niobate: the role of the 

delayed Raman response 
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We review measurements in the literature regarding Kerr nonlinearity of lithium niobate. Of 
particular interest is to separate the value of the electronic Kerr nonlinearity and the strength of the 
delayed Raman nonlinearity. The spectral shape of the Raman response in LN is discussed and we 
restrict ourselves to the case where the pump is polarized along the optical axis. This configuration 
is important for so-called type noncritical interaction, which is used in QPM devices and in recent 
supercontinuum generation and few-cycle soliton compression experiments. 



I. INTRODUCTION 

Lithium niobate (LiNb03, LN) is one of the earliest 
frequency conversion crystals that was used, and its op- 
tical linear and nonlinear properties have been studied in 
a large number of publications over the past 5 decades. 
It remains today one of the most important near-IR non- 
linear crystals, so it is striking that its cubic nonlinear 
properties still need to be revised and investigated, as 
they play a crucial role in recent ultrafast experiments 
[1—3]. In fact, despite very early and thorough studies 
of the incredibly complex and rich nature of the cubic 
nonlinearity, where the Raman spectra were measured in 
detail [4-6] , still many key parameters needed for model- 
ing and understanding ultrafast interaction are unknown. 
The purpose of this review is to present various measure- 
ments of the Kerr nonlinearity in LN and through the 
analysis try to propose good parameters usable of mod- 
cling the instantaneous (electronic) and delayed (Raman) 
nonlinearity in LN. We focus on the specific case of so- 
called type interaction, where both the pump and the 
generated second-harmonic are polarized along the opti- 
cal axis; this interaction is widely used in the nonlinear 
optical community as it exploits the largest d.33 compo- 
nent of the quadratic nonlinearity, and it is in particular 
used for devices based on quasi-phasc matching (QPM) 
interaction. 

An initial motivation for carrying out a literature study 
was a recent experiment by one of us: In Schiek et al. [7] 
four-wave mixing (FWM) was used to measure the Ra- 
man response of LN. A strong pump at X p = 1.064 fxm 
and a weaker signal varied between A s = 1.07 — 1.25 /im 
were overlapped in the LN sample. This generated a 
new photon through FWM with the frequency 2w p — lo s . 
The LN was X-cut, and the pump and probe were both 
polarized along the Z-axis, which contains the largest 
quadratic nonlinear tensor component (I33. This corre- 
sponds to 9 — 7r/2, so the beams have extraordinary po- 



larization. By analyzing the FWM signal, it could be bro- 



ken down to = Xb^ " 
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Xm- Thus, contri- 



butions to the total nonlinearity came from the electronic 
Kerr effect Xe? , cascaded quadratic nonlinearities Xcasc 5 

(3) (3) 

the Raman effect X r ■, an d finally IR cascading x\r- By 
using well-known formulas for cascaded quadratic nonlin- 
earities, tabulated Raman data from the literature, and 
IR dispersion relations, the electronic component could 
be extracted. The surprising result of this study was that 
it found a very large electronic Kerr nonlinearity for the 
ZZ ZZ component 



X [ e = 0.56 x 10 4 pm 2 /V 2 
n\ E = 34.6 x 10~ 20 m 2 /W 



(1) 
(2) 



When modeling femtosecond pulse propagation the "to- 
tal" Kerr nonlinearity is usually given as the sum of elec- 
tronic and the Raman term, and in this case it gave 



= 1.4 x 10 4 pm 2 /V 2 



2, tot 



83.3 x 10~ 20 m 2 /W 



(3) 
(4) 



and we can therefore calculate the overall Raman fraction 

as 



fn = Xr 3 Vxt ( ot = 0-58 



(5) 



which is a very important number to know when im- 
plementing the Raman effect in a numerical simulation. 
Finally the following ratio was found 



Xr ' /max\xR (peak)} 



=0.165 



(6) 
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This represents the ratio between the DC value and the 
peak of all the Raman lines (i.e. the value at zero fre- 
quency when the Raman spectrum is normalized to have 
its peak at unity). 

These numbers are much larger than what the litera- 
ture have found so far: typically Z-scan measurement re- 
vealed near- IR values of n\ ~ 10 x 10~ 20 m 2 /W for this 
tensor component [8] , and although it was not specifically 
stated there this value is the sum of electronic and Raman 
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nonlinear contributions (sec also review in [9]). Moreover, 
consensus from the fiber optics community says that Ra- 
man fractions should be on the order fn = 0.1 — 0.2. The 
question is whether LN literature data can confirm the 
surprisingly large values found by Schiek et al. 



II. RAMAN TERM 

In the standard definition of a vibrational Raman- 
mode, it is described in time domain as a delayed frac- 
tional response adding to the instantaneous electronic 
cubic nonlinearity. Consider the nonlinear susceptibility 

[10] 



X W (t-tl,t-*2,t-*3) 



xif t R(t-h)6(t-t 2 )6(t-t 3 ) (7) 



where Xtot i s the "total" nonlinearity; as we see below it 
is the sum of the electronic and DC Raman nonlinearity. 
The response function R is defined as to include contribu- 
tions from both electronic origin and from the vibrational 
Raman effect from a number of Raman modes 



R(t) - (i - f R )S(t) + f R h R (t) 



(8) 



where fa is the fractional contribution of the Raman re- 
sponse to the total nonlinear polarization, and 



M*) 



(9) 



where fn. a is the relative strength of mode a so 
So- fn,v = 1- -R is normalized so J_ oo dtR(t) = 1, and 
thus hn is normalized as well. 

In frequency domain we have equivalcntly 



(3) 



where 



(3) 



(10) 



(11) 



and /ii? jCT (fi) is the Fourier transform of h,R ylJ {t). This 
definition means that the DC value of the Kerr nonlinear 
coefficient is 



(12) 



This is the value most often measured, since standard Z- 
scan measurements have been performed with multi-ps 
pulses (see overview in Ref. [9]). 

The standard model for h R is a complex Lorentzian in 
frequency domain 



l R,a 



(13) 



whose normalized time form is [11] 

T 2 + T 2 

hR,a(t) = — 2-^- ex P(- t / T 2.<r)sin(i/Tl, C r), t > 

(14) 

where ri )(r is the period of the temporal oscillations of 
mode a and r-x,a is characterizes its decay. Note that the 
response function is only defined for t > for reasons 
of causality. Taking the Fourier transform we get in fre- 
quency domain 



/oo 
dte- im h R , a (t) 
-oo 



'l.cr T- '2,er 



(15) 



The normalization is checked by noting that hn^i^l = 
0) = 1 and h R)(T (n = 0) = J^dthR^t). Now let us 
rewrite the part of the Raman response of Eq. (13) that 
contributes to the frequency variation, which we conve- 
niently normalize by considering that at 51 = it should 
be unity 



O 2 4- O 2 



'1,<T ' '2,<T 



(16) 
(17) 



where in the first line we have introduced 

ft 2 , CT = r ff /2 (18) 
n 1>a = yffi viz a = y^ 2 ni a (19) 

and now fl2,a is the HWHM linewidth. In the second line 
we used 



n,„ = n£ = [ l 4-il/A]- 1 * 

T2,a = ^2,ct = 2/r CT 

and we now have the desired form of Eq. (13). 



(20) 
(21) 



III. RAMAN PARAMETERS AND RAMAN 
GAIN MEASUREMENTS OF JOHNSTON AND 
KAMINOW 

Early studies by Johnston and Kaminow [5] found four 
Raman lines, and even measured an absolute scattering 
strength for each mode. Using these data we can recon- 
struct the absolute values of the Raman peaks in mks 
units using the equations of [6]. The only uncertainty 
is that they are unclear about the reported linewidths 
(whether what they report is the FWHM or not), and 
compared to other studies in particular [5] seems to mea- 
sure quite narrow linewidths, whereas [6] only measured 
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the two dominant lines (at different temperatures) with 
results more in the trend of other studies. In Table I of [ ] 
the transverse "T" modes of the type interaction (all 
waves polarized along Z) are reported with scattering ef- 
ficiencies S33/MO in units cm sr . The linewidths are 
reported as 2~f a , which we assume to be FWHM (i.e. 7 
is the HWHM, and the notation 7 is used to emphasize 
that it is different from the T used above: 7 = T/2), how- 
ever the linewidths of 2-f a = 15, 8, 6 and 15 cm -1 for the 
four modes located at 253, 275, 334 and 637 cm" 1 , re- 
spectively, seem quite low. Not unsurprisingly, they later 
correct these linewidths, stating that they were 30% too 
small, see [12, p. 3497] and [6, p. 1049]. Proceeding with 
specific calculations, Eq. (7) in [6] can be used to find the 
power-gain coefficient g s relative to the pump intensity 
I p (note the corrected equation from the erratum [13] is 
used) 



o 




9s 

In 



16tt 2 c 2 



S a /ldn 



hujj n 2 (n 0)<7 + 1)27„ 



(22) 



where (all in mks units) lo s is the absolute frequency 
of the Stokes wave, n p and n s are the linear refrac- 
tive indices of the pump and Stokes waves, no t(T = 
l/[exp(fiw -/fc^r) — 1] is the Bose population number 
of the Raman frequency at the temperature T. Finally, 
Scr/ldCt is the scattering efficiency (again the corrected 
equation from the erratum [13] is used) 
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\(da/dq) Q \ 



(23) 



where p a is an effective reduced mass per unit volume 
of the mode. An important parameter (da/dq)o is intro- 
duced, which is the derivative of the optical polarizability 
with respect to the normal coordinate q of the vibration. 
Below we discuss this in further detail and show how it 
is connected with the first order transition hyperpolariz- 
ability. Note also the well-known scaling with . 

The experiment was conducted at room temperature 
and with X p = 0.488 /im. We now calculate the peak gain where by definition n 
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FIG. 1. Data from Kaminow and Johnston [5] measured for 
A = 0.488 fim. The Raman response Xr vs. inverse wavenum- 
ber (top) and the corresponding Raman contribution to the 
nonlinear refractive index n| r vs. angular frequency (bot- 
tom). The blue curve is the real (dispersive) part while the 
red curve is the imaginary (absorptive) part. The dashed black 
line indicates the absolute value. Note that the linewidths 
have been increased 30% to the values of [5], see more in [12, 
p. 3497] and [6, p. 1049], and that the gain calculations use 
the corrected equations from the erratum [13]. 



values of the four modes, which gives [5, 13] 

Y = {14.7, 7.2, 1.7, 13.5} cm/GW (24) 
Ip 

(using 30% larger 7^) 

Now we can relate this Raman gain to a nonlinear sus- 
ceptibility by using that the polarization response of a 
Kerr interaction in the type case is [Eq. (B8) in [!)]] 



no + n^I, i.e. a nonlinear index 
change that is proportional to the intensity. From the 
SVEA we get that the cw wave equation for the Stokes 
field £ s at the frequency uj s = uj p — uj a is 



2ik, 



dE s 
dz 



— 2^NL 



e c 



(27) 



P, 



(3) 



NL 



Si 



(25) 



where k s = n s oj s /c. The Stokes wave sees Raman gain 
from the XPM of the pump, and under the approxima- 
tion of a weak gain the Stokes SPM is neglected. Thus, 
neglecting SPM 



(3) 

where X c g a is the effective nonlinear tensor component 
that depends on the birefringent properties of the crystal; 

for type we have Xctt a = Xe&,ij = ■ The connection 
to the Kerr nonlinear index is 



d£ s 
dz 



. 3w s 
t 

An s c 



^\uj a )\£ p \ 2 £ s 



(28) 



3x (3) 
4en en 2 



(26) 



The gain g s in [6] is defined as the power exponen- 
tial gain e 9sZ , so we look for solutions = g s I s , 
where I s = eon s c\£ s \ 2 /2 is the intensity. Since H^- = 



4 



he n s c2Re[£, 



d£ 3 
a dz 



we get that 
3u? 



2n s c 
3u> 



n s n p E c 
(in 



Im[ X ( 3 )(^)]|f p | 2 
hn\x% (w s )]I p 
Im[XR {u s )]Ip 



\ s n s n p e c 



(29) 
(30) 
(31) 



where we define the Kerr nonlinearity in frequency do- 
main from Eq. (10) as x (3) H = Xe? + Xr\u)- Eq. (30) 
is identical to [14, Eq. (7.36)]. Since the Raman suscepti- 
bility is imaginary and negative at the Raman resonance 
frequencies, the Stokes wave experiences Raman gain, as 
expected. Also note that since the gain is given by the 
imaginary part of the nonlinear susceptibility, then the 

(3) 

electronic part \e wm n °t affect the result as it is real. 
We can now use this result to calculate the imaginary 
Raman peak values as 



Xfl,l(peak) = |Im[x^K = 
_ \ s n s n p Eoc g s 



(3), 



Gtt 



(32) 



Eq. (32) is a factor of 2 smaller than Eq. (10) in [2]; 
it could be because they calculate the amplitude gain 
= a s £ s [this is at least the approach used in [15], cf. 
Eq. (10.3.21) there - but note that he uses a completely 
different definition of the electric fields and susceptibili- 
ties], and as is well known a s = g s /2. Carrying out the 
calculations, we obtain the peak values [5, 13] 

xgl(peak) ={5.2,2.5,0.61,4.8} x 10 4 pm 2 /V 2 , (33) 
(using 30% larger j a ) 

Assuming again a Lorentzian-based response for the Ra- 
man we can write it on the form 



y OW _ (3) "a 



cr=l 

4 

Xr ] f R <°~, 

<T = 1 



i2jM - n 2 



(34) 



(35) 



where as usual \r ' = S CT Xr ' a an< -f wc navG used the 
requirement — Im[x^ (w ff )] = \^r \ (peak) to define 



(3) 



X^l(peak) 



fR,a 



(3) 
XRa 

(3) 
Xr 

,( 3 ) 



UvUpljXRj 



a (^-^y- + (2^ni 



(36) 
(37) 



Note that knowing Xr CT (P ea k) we can find the individual 

(3) 

Xr u by solving the system of equations for all a that 



Eq. (36) defines. These definitions ensure that the peak 

values match the calculated ones Xfl^peak), and that 
everything is normalized as before. The calculations give 
[5, 13] 



X ( r = 0.62 x 10 4 pm 2 /V 2 

n 2,R = 34.3 x 10~ 20 m 2 /W 

f Ra = {0.635,0.105,0.020,0.240} 



(38) 
(39) 
(40) 



Note that these values are quite insensible to the 
lincwidths used. The time constants are [5, 13] 



ri, CT = {21.0,19.3,15.9,8.3} fs 
r 2)CT = {544, 1021, 1361, 544} fs 
(using 30% larger 7 CT ) 

Finally we mention the ratio [5] 



(3) 
Xr' 



xk 3 l(P eak )] 



=0.119 



(using 30% larger -f a ) 



(41) 
(42) 



(43) 



This represents the ratio between the DC value and the 
peak of all the Raman lines (i.e. the value at zero fre- 
quency when the Raman spectrum is normalized to have 
its peak at unity). 

Let us address the question about the measurements 
performed in the subsequent paper [6, 13], where the 256 
and 637 cm -1 modes at 300 K were reported to have 
linewidths of 2j = 23 and 20 cm -1 , respectively, with 
the same reported scattering efficiencies. Inserting the 
numbers of Table II of [6] in Eq. (22) we get [6, 13] 



g s /I p ={12.5, 12.9} cm/GW 



(44) 



at 256 and 637 cm -1 , respectively, in agreement with 
Table II of [6] when correcting for the factor of 2 mistake 
[13]. Using the same approach as above the peak values 
of the two modes are calculated to be [6, 13] 



xSl(peak) ={4.4,4.6} x 10 4 pm 2 /V 2 



(45) 



These values agree quite nicely with their previous re- 
sults, Eq. (33), when all the corrections are applied. We 
can still use the value from Eq. (38), as calculated from 
their previous measurements, as the DC value, because 
only the linewidths separate the measurements, and as 
mentioned above they do not influence the DC value. 
Assuming the stronger of these two lines to be the peak 
of all the Raman lines, the ratio of the DC value and the 
peak of all the Raman lines is now [(>, 13] 



xk 3) /max[x^ (peak)] =0.133 



(46) 



again close to the previous result Eq. (43). 

Johnston [16] measured the gain at A = 1.064 /xm of 
the v a = {256,637} cm" 1 lines as 



g s /I p = {5.75,6.0} cm/GW, 
(correcting for a factor 2) 



(47) 
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Calculating the peak values of the Raman susceptibility 
and assuming the same linewidths as in [G] (although it 
is not specifically defined) we get [Hi] 



Xfl.l(peak) = {4.1,4.5} x fO 4 pm 2 /V 2 
(correcting for a factor 2) 



(48) 



IV. 



OTHER RAMAN MODELS AND RAMAN 
GAIN MEASUREMENTS 



A measurement of the two main Raman gain peaks 
is reported by Boyd [15, p. 480], where the v a = 
{256,637} cm" 1 lines are reported to have linewidths 
Ai/ a = {23,20} cm , and peak gain 



g s /I p = {8.9,9.4} cm/GW, 



(49) 



measured at 694 nm. Calculating the peak values of the 
Raman susceptibility and assuming that Av a = T a (al- 
though it is not specifically defined) we get 



Xfl,l(peak) - {4.2,4.6} x 10 4 pm 2 /V 2 



(50) 



Finally, in Chunaev et al. [17] measured the gain at 
A = 1.047 /im of the 256 cm" 1 mode 



g s /I p = 5.1 cm/GW 



The decay time of the mode was found as t-i = 0.38 ps, 
which gives a lincwidth of T = 28 cm -1 . Calculating the 
peak values of the Raman susceptibility we get 



Xi?,l(P eak ) = 3 - 6 x 104 pm 2 /V 2 



(52) 



We include for completeness the Raman spectrum of 
Phillips et al. [2] in Fig. 2. They found 7 modes, the first 
four are the standard ones plus 3 minor modes. They 
reported a DC ratio 



(3) 

Xr 



max [x«,l(P eak )] 



0.17 
0.156 



(53) 
(54) 



where the latter value was found from calculating the ra- 
tio from their Table 1. The 3 extra modes they include 
do not contribute much to this ratio. The chose to use 
maxfx^ (peak)] = 5.3 x 10" 21 m 2 /V 2 for A = 1.043 ^m. 
This was determined based on data from numerical sim- 
ulations. We then get using the DC ratio 0.17 



(3) 

Xr 

n 2,R 



0.090 x 10 4 pm 2 /V 2 
5.5 x 10~ 20 m 2 /W 



(55) 
(56) 



When comparing to the previous measurements these 
values are very low. Note also how strong the large- 
frequency mode is relative to the main low-frequency 





0.4 








0.2 


a 




ex 

-t- 


0.0 


o 






-0.2 


So; 






-0.4 









- / 




A \ 


l : 
\ : 


■J 






\ 




1 







-500 



500 



10 

-10 
-20 
-30 



y[] 


L/cm] 


u 




1 


V 


i 

1 : 
\ ; 
\ 








Y ~\ 





-200 -100 100 200 

fl [rad/ps] 



(51) 

FIG. 2. Data from Phillips et al. [2]. The Raman spectrum 



is scaled to give the peak value of max[xfl' CT (peak)] = 5.3 x 
10" 21 m 2 /V 2 of the largest frequency mode, applicable for 
A = 1.043 /urn. 



mode. Finally, note that the spectral form they use 
a.j/(fj — f 2 + 2ijjf) is not normalized. 

Dclfyett [18] is again another story, because they fit to 
what they claim to be a Lorentzian and get the Raman 
contribution to the total nonlinearity of 



H( U ) = £ a - 



(57) 



This form is problematic because (a) it is not (conjugate) 
symmetric towards negative frequencies (b) the real part 
is negative at zero frequency. We remark that one can 
reduce the classical Lorentzian form Eq. (13) to this form 
close to a resonance. Let us transform their form into the 
standard form, Eq. (13) by normalizing it. Thus, with the 
parameters 



do 



{0.61,0.27,0.12} 



(58) 
(59) 



and obviously 2"j a = T a we can use the standard form. 
We have in Fig. 3 used as a a and 7 CT the parameters a, 
and r, in Table I of [181. We also find the DC ratio to 
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FIG. 3. Raman spectrum from Delfyett et al. [18] measured 
for A = 0.532 fim, and using a; as the oscillator strengths. 
Top: their model, Eq. (57). Bottom: equivalent model using 
the symmetric Lorentzian Eq. (13). Notice that the real part 
of the top plot is negative in the model they use. 



peak as 



(3) 

Xr 



max [Xflt(peak)] 



0.075 



(60) 



V. FREQUENCY DEPENDENCE OF THE 
ELECTRONIC AND RAMAN 
SUSCEPTIBILITIES 
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FIG. 4. As Fig. 3 but using a/\n) as the oscillator strengths. 
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FIG. 5. Left: The degenerate G 2 parameter [20, p. 279] of 
LN vs. wavelength, calculated using E g — 3.9 eV [8]. To a 
good approximation the frequency dependence of the elec- 
tronic nonlinearity n^g oc G2. The vertical line at 636 nm 
marks the wavelength of a photon carrying half the bandgap 
energy. Right: calculated n\ E from Eq. (61) using the G2 in 
the top graph. 



Now we have Xr > but in order to estimate the 
Raman fraction fn we need to know the total Kerr 
nonlinearity (or the electronic Kerr nonlinearity), and 

fn = XflVxtot- At this point it is relevant to dis- 
cuss the frequency dependence of the nonlincaritics. 
We know that the nonlinear susceptibilities change 
with frequency. Miller [19] first pointed out that 
5 = x (2) ( w i ; w 2,w 3 )/[x (1) ( w i)x (1) ( w 2)x (1) ( w 3)] is roughly 
constant, hence it is known as Miller's 5. One can there- 
fore given the value at one wavelength apply Miller's scal- 
ing to estimate the nonlinear coefficient at a new wave- 
length. Such a relation also holds for the higher-order 



nonlinearities. It can be derived from treating the atom 
and the electron using the Lorentz model of a classical an- 
harmonic oscillator, and then including a nonlinear force 
of the electron (see e.g. [15]). Therefore one can only use 
this approach to scale the electronic nonlinear suscepti- 
bility. 

A more accurate model is the two-band model, which 
is [20, Eq. (55)] 

n{ E (u;)(SI) = K' G 2 (tkj/E g ) (61) 
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where K' is a material constant, fitted to K' = 7.3 x 
10~ 9 eV 35 m 2 /W in wide-gap dielectrics [8], and E p = 21 
cV is the Kane energy, which is constant for most mate- 
rials. With these two constants one can a priori estimate 
the electronic nonlhiearity, in particular the enhancement 
close to the two-photon absorption regime. The curves 
are plotted in Fig. 5. 

The question is now: does the DC Raman nonlinearity 

(3) 

Xft scale with wavelength? And if yes, how does it scale? 

In order to answer this, consider a classical description 
where Raman scattering can be modeled as a simple har- 
monic oscillator. Consider the standard case (no Raman): 
the electric field induces a change in the molecular dipole 
moment p — e^aE by shifting the electrons away from 
the nucleus (consistent with before we use tilde to denote 
the choice where Eq is kept apart from the polarizability, 
as we below have another case where this is not done). 
The linear susceptibility is given by x = Na(u), where 
N is the atom number density, and the polarizability has 
a Lorcntz form 



a(ui) 



1 



e 2 /m 



e w 2 . - u 2 - 2i^ a u: 



(62) 



where e is the electron charge, m the electron mass, uj a 
the dipole resonance frequency and "f a the dipole damp- 
ing rate. This approach assumes that a = ao is constant 
in time, and therefore the electrical-field shifted electron 
cloud is in equilibrium at its new shifted position away 
from the nucleus. The Raman effect occurs when the po- 
larizability has a temporal dependence around the equi- 
librium value [15] ao as a(t) = uq + q(t)(da/dq)o, where 
(da/dq)o is the (first) hyperpolarizability [21, Ch. 15], 
also called the differential polarizability [14, p. 199], and 
where the subscript '0' indicates that it is a static value at 
the particular applied optical frequency. The Stokes field 
turns out to have the following nonlinear susceptibility 
[15, Eq. (10.3.19a)] 



xk 3) K) 



e (N/6m)(da/dq)l 



cj 2 - {ui p - uj s ) 2 + 2i7 (T (w p - uj s ) 



(3) 

Thus, at zero frequency we have Xr K = w p ) = 
£f)N{da/dq)^/{QrnLLi 2 ). The only term that could have 
a frequency dependence is the (da/dq)o term, but in the 
classical model it is a purely phenomenological quantity 
so we cannot from this approach say anything. Certainly 
a depends on frequency, this is stated in Eq. (62). But is 
the hyperpolarizability (da/dq)o dispersive? According 
to [14, p. 199] the Raman susceptibility close to a res- 
onance can according to a classical derivation be stated 



N\(da/dX) \ 



6feo(w CT - uj p + lo s + i-ja-) 



(64) 



where (da/dX) is the differential polarizability, and as 
the notation indicates it is defined in a different way than 



(da/dq)o, see more below. In the same limit Eq. (63) 
becomes [15, Eq. (10.3.19b)] 



e Q N(da/dq) 2 



12mw cr (w - - Up + uj s + i"f a ) 



(65) 



Care must be taken in comparing these two results: in 
[14] the dipole moment p induced by an electric field E 
defines the polarizability as p = aE. Instead in [15] it is 
defined as p = eoaE; this ensures that a has units [m 3 ] 
and that x*- 1 ) = JVa(u) holds. We of course have simply 
that a = Eqol, and that [a] = Fm 2 . More importantly, in 
[14] A is a dimensionless coordinate, so (da/dX)o has 
the same units as a. This means that the (reduced) elec- 
tron mass m used above is replaced by S/2w CT , stated in 
[14] as a reduced mass but really is a reduced mass nor- 
malized to a characteristic length suitable for normal- 
ization of A. In essence if one makes the substitution 
y // h/2muj rT X = q the two results become identical. In 
[14] the classical derivation leads to the following identity 
(da/dX)o = aj g (—uj s ;ujp), linking the result to a quan- 
tum mechanical derivation where the first order transi- 
tion hyperpolarizability, describing the generation of a 
Stokes wave, was introduced as [14, Eq. (4.115)] 



i%(-u e ;u p ) = ^12 



r fi e P ' r gi 



n 



(66) 



where the subscript fg symbolizes the the transition from 
a molecular ground state g to a excited final state / with 
transition frequency f2/ ff , and the sum over the index i 
refers to all the final states as observed from the molec- 
ular ground state g. We have here used the superscript 
'tr' to emphasize that this transition hyperpolarizabil- 
ity is different than the molecular polarizabilities a (or 
a). They are defined in similar ways but with different 
scopes: a is the first order expansion coefficient in the 
dipole moment when Taylor expanded around the applied 
electric field (e.g., p = p + aE) [14, Eq. (4.83)], while 
(63) ^tr j g |- ne £ rs -j. or( j er expansion coefficient in the polar- 
izability when Taylor expanded around the vibrational 
position coordinate (q or A) from its equilibrium value 



(e 



g., a = a 



a 



r A) [14, Eq. (6.121)]. The electric field 



is defined as = ejEj, so is the dimensionless unit 
vector. The vibrational coordinates of the Raman modes 
are given by Vj , and e is the charge of the electron. It is 
important to note that [aj g (— w s ; uj p )] = Cm 2 /V = Fm 2 , 
i.e. the same as a and (da/dX) . We can now compare 
the two cases and get 



(da/dq) Q 



;w P )| 



(67) 



and obviously (dot/dq) a = Eo (da/dq) Q . This shows that 
we can expect the differential polarizability to experience 
resonant enhancement when the pump photon lo v or the 
Stokes photon lo s comes close to a molecular resonance 
frequency. [Eq. (66) is derived assuming interaction far 
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from such a resonance, and thereby the damping terms 
have been removed from the denominators in the square 
brackets in Eq. (66).] The numerators should actually 
express quantum mechanical expectation averages over 
the ground and final stage wave functions and the elec- 
tric dipole moment operator. Whether these quantities, 
which could be calculated in the simplest case using e.g. 
variation-perturbation theory, vary with pump frequency 
is an open question. 

We can now try to express the peak gain when ui p — 
lu s = ui a using the polarization derivative. The peak Ra- 
man susceptibility Eq. (65) becomes at this frequency 



Xr (peak) 



Z N{da/dq)l 



(68) 



We can now use Eq. (32) to get the peak Raman gain 
coefficient as 



gs_ _ uj s N(da/dq)l 
I p Amuj^nsnpC 2 



(69) 



If we do the same by inserting Eq. (23) into Eq. (22) 
then we get the same result, provided that we substitute 
(da/dq) = e (da/dq) . 

Hcllwarth used the Born-Oppenhcimcr approximation 
to derive the nuclear (vibrational) contribution to x*- 3 ' 

and argued that it must scale as x^r X < " 1 H UJ p)x^ 1 H L °s) 
[22, Sec. 13.2.1], which then implies 



X^oc^- 1)^-1)^(^-1^ 



(70) 



which means that the Raman gain scales as 

Im tx« ]/(VvO, ci - Ec f ( 31 )' or 



g s /I p oc (n^ - l){n 2 s - l)/(\ s n s n p ) 



(71) 



This presumably is what leads to the simple scaling- 
factor of the Raman gain (n 2 — l) 2 /n 2 reported several 
places [23-25]; what Stolen [23] refers to when quoting 
the (n 2 — l) 2 /n 2 scaling is the gain scaling besides the 
inverse scaling with the wavelength and he then presum- 
ably simplifies n p ~ n s . However, it also seems a sort of 
conclusion from [24, 25] that this scaling alone does not 
explain the measured increase in the "normalized" Ra- 
man gain (normalized to scale out the A" 1 dependence) 
when the pump frequency is close to an electronic tran- 
sition (i.e. the UV loss edge), see e.g. [24, Fig. 4]. The 
increase was instead explained as a resonance enhance- 
ment [24, 25]. 

How can we model this resonance enhancement of the 
Raman susceptibility? Eq. (66) shows a resonance en- 
hancement, but this is connected to a proximity to the 
molecular transition frequencies, and not electronic tran- 
sition frequencies. Levenson and Bloembergen [26] ex- 
pressed the Raman terms with the local-field corrected 
polarizability a R , which corresponds to the first-order 
transition hypcrpolarizability as defined in Eq. (66); this 
can readily be checked since [26, Eq. (8)] corresponds to 



Eq. (67) when neglecting what presumably is the local- 
field correction term [(n 2 + 1)/3] 2 (although the standard 
field-correction factor is [(n 2 + 2)/3] 2 , see [27, Eq. (4)]). 
Moving on, they introduced a phenomenological scaling 
as [26, Eq. (18)] a R cx 1/(£L - tjw p ), where £L is the 
is the energy of the direct band gap. They do mention 
that the resonance energy is chosen rather arbitrarily, 
and could also be chosen to coincide with that of the 
Sellmeier equation. For LN we have E g = 3.9 eV, corre- 
sponding to A = 0.32 fj,m. The UV pole of the Sellmeier 
equation for n e instead lies around A = 0.20 — 0.21 /xm. 
Kato and Takuma [27] relied on a different quantum me- 
chanical derivation of the Raman response by Peticolas 
et al. [28], and they arrived at a resonance enhancement 
of what corresponds to the transition hypcrpolarizability 
as R oc [(uj 2 — uj 2 )(uj 2 — uj 2 )]' 1 , where wo is the energy 
difference between a ground and an excited electronic 
state without coupling to molecular vibration [27, Eq. 
(3)]. The total Raman gain was then given as oc |i?| 2 . 
This suggests an enhancement of the Raman gain that 
scales as [(uj 2 — ui 2 )(ui 2 — a; 2 )]" 2 , assuming a symmetric 
excitation of the vibration (I assume this means stimu- 
lated Raman scattering, i.e. that the two pump photons 
excite the same Raman mode, and thus the degenerate 
case of [14, Fig. 4. 5(b)]). Other authors have without pro- 
viding any reasoning chosen to model the scaling of the 
Raman gain coefficient as [29-31] 



il =D 

I P (y% 



v 2 ) 2 
p> 



(72) 



where vq is the resonance frequency. Similar results were 
also found in [32, Eq. (2)]. In fact, there they find using 
the Born-Oppenhcimer approximation that the transi- 
tion hyperpolarizability a = A + B, where the param- 
eters A oc (v 2 + v 2 )/(v 2 — v 2 ) 2 , corresponding to a diag- 
onal "symmetric" excitation, and B oc (v e v s + v 2 )/(v 2 — 
v 2 )(v 2 — v 2 ), corresponding to an off-diagonal "asymmet- 
ric" excitation, and v e and v s are virtual electronic states 
(i.e. energy differences between ground and excited elec- 
tronic states). The total Raman cross section then scales 
as iA*|a tr | 2 so they introduce two dimcnsionlcss frequency 
factors 

F A = vl(vl + vl)/{vl-vlY (73) 
F B = 2v 2 s (v & v s + V l)l(v 2 e - V 2 v )(v 2 s - vl) (74) 

Depending on the nature of the experiment, the Raman 
scattering cross section can now scale as F\, correspond- 
ing to a symmetric excitation, as F^, corresponding to 
an asymmetric excitation, or a mixture, (Fa + Fb) 2 ■ 

In Fig. 6 we plot the gain data vs. wavelength, and we 
found surprisingly that the resonant enhancement scaling 
laws did not fit very well with the data. Concerning the 
Raman gain, left plot, we found that the simple scaling 
found by Hellwarth, Eq. (71), seems to fit. However, a 
large part of the wavelength dependence is attributed to 
the 1/A S dependence, so a more accurate way of under- 
standing the frequency dependence is to plot the peak 
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FIG. 6. Top left: Raman Stokes gain g a /I p as measured at 
various wavelengths. The red (black) dots are the 256 cm" 
(637 cm -1 ) mode, and the curves follow the scaling suggested 
by Eq. (71), and is chosen to be centered around the mea- 
surement at 694 nm. Top right: the same measurements ex- 
pressed through the peak Raman susceptibility \^{peak). 
The curves now follow the scaling suggested by Eq. (70). Any 
resonance scaling was not found to match the data. Bottom: 
suggested frequency independent parameter Sr, cf. Eq. (75) 
vs. wavelength. The error bars are indicating an estimated 
10% standard deviation. 



VI. MEASUREMENTS OF THE KERR 
NONLINEARITY 

In DcSalvo et al. [8] they used a Z-sc&n. technique to 
measure nf = 48 ± 7 x lCT 14 csu (n| = 9.1 ± 1.3 x 
10 -20 m 2 /W) at Ai = 1.064 /xm in an A-cut congru- 
ent LN crystal (undoped) with the beam polarized along 

(3) 

Z, which means they measured the Xzzzz tensor com- 
ponent. They assumed that no cascading contributed to 
this value, since the phase mismatch is large. We cal- 
culate using the Sellmeier equation for congruent LN of 
[36] that rie(cJi) = 2.1558 and n e (2u 1 ) = 2.23421, so 
Afc = 926 mm -1 . We now use the well-known quadratic 
cascading Kcrr-like nonlinear refractive index [37] 



\pE0cnfn2Ak 



(76) 



where d c g is the effective quadratic nonlinearity in units 
[m 2 /V 2 ] and Afc = fe— 2fci, kj = uijUj/c, and rij, j = 1,2 
are the linear refractive indices of the pump (j = 1) and 
the second harmonic (j — 2) With d c s = ^33 = 25.2 
pm/V [38] we get from Eq. (76) that n£ ca8c = -29.4 ± 
2.9 x 10~ 20 m 2 /W; the uncertainty stems from a ±5% 
accuracy in d 33 [38, p. 2272]. Thus, cascading contributes 
quite strongly to the measured nonlinear refractive index. 
A 30 ps pulse duration was used, so they measured the 
electronic plus the DC value of the Raman term, see also 
Eq. (12). If we subtract the estimated cascading value, 
the actual total Kerr nonlinear index would be 



susceptibility, which we have done in the right plot. The 
red data, i.e. the low-frequency Raman mode, seems to 
scale like suggested, i.e. as Eq. (70), while the black data 
from the high-frequency mode seems not to depend on 
wavelength at all. In order to remove any wavelength de- 
pendence one could suggest a sort of Raman- version of 
the Miller's delta as 



1 = 38.6 ±3.2 x 10" 20 



Sr = 



(3) 

Xr 



X«H) X «(a, s ) 



(75) 



which then should remain constant with wavelength. We 
have plotted these data in the bottom figure, and obvi- 
ously the scaling seems reasonable for the low-frequency 
mode, while the high-frequency mode seems to increase 
with wavelength. It should here also be noted that for 
some strange reason there is no agreement in the litera- 
ture whether the low-frequency mode should be stronger 
than the high-frequency mode [4, 5, 7, 33-35], or vice 
versa [2, 6, 16, 18]; note in this connection that above 
room temperature the high-frequency mode becomes 
stronger [35]. This uncertainty in determining the mode 
strengths could influence the analysis of the data. Again 
we stress that we did not find any resonance enhance- 
ment for the data, which could be because the possible 
resonances lie too far into the UV to be observed. 



n 2,tot 



m 2 /W 



(77) 



where the error has been calculated using propagation of 
uncertainty cr t 2 ot = a%_ scan + <J^ SC ; the relative error is 
8.4%. 

DeSalvo et al. also performed the Z-scan measurement 
at A = 532 nm and measured nf = 440 ± 70 x 10~ 14 esu 
(n| = 83 ± 13 x 10~ 20 m 2 /W). The cascading contribu- 
tion at this wavelength is small, but not insignificant. 
To estimate it we calculate using the Sellmeier equa- 
tion for congruent LN of [36] that n e {u)i) = 2.23421 
and n e (2wi) = 2.28590, so Afc = 14,759 mm" 1 . We 
then use Miller's scaling to translate d c s = 1^33 = —25.7 
pm/V measured at Ai = 0.852 /xm [38] to 532 nm and 
get d e ff = —49.6 pm/V. Thus using Eq. (76) we have 
n Lasc = -10.4 ± 3.1 x 10~ 20 m 2 /W. The uncertainty of 
this number lies mostly in the d e g value (we estimate at 
least 10% error from Miller's scaling on top of the 5% ac- 
curacy of the measurement itself), because although the 
Sellmeier linear refractive index at 266 nm is not guaran- 
teed to be that accurate, we checked that UV measure- 
ments [39] support the calculated Sellmeier value. The 
total Kerr nonlinearity is then when correcting for cas- 
cading 



93 ± 14 x 10" 



m 2 /W 



(78) 



Such a high value seems surprising, but consider that 
the wavelength corresponds to a photon energy of 2.33 



10 



200 
100 



-100 

-200 



1000 2000 3000 




1000 2000 3000 





4 












1 









v[THz] 



-1000 1000 2000 3000 
t[fs] 



-, 200 
r L. 100 

° o 




4 








a -100 
= -200 
















v[THz] 



-1000 1000 2000 3000 
t[ft] 



FIG. 7. Convolving the Raman response of Fig. 1 with a Gaus- 
sian pulse having intensity FWHM of 1000 fs (top). 200 fs 
(2. plot), 100 fs (3. plot) and 50 fs (bottom). The left plot 
shows the Raman response, and the Fourier transform of the 
Gaussian intensity (displayed using arbitrary scaling in the y- 
direction) is shown as a dashed line. The right plot shows the 
convolved time response (real part as blue and imaginary part 
as red, the latter always being zero as the pulse is symmetric 
around zero frequency) , and the black dashed line corresponds 
to the long-pulse limit, where n^n(t) — ► n^l^t)! 2 . 



fore concluded that the measured nonlinearity is purely 
instantaneous. We checked this by convolving the Raman 
response of Fig. 1 with a Gaussian pulse having intensity 
FWHM from 50-1000 fs, see Fig. 7. What we are inter- 
ested in is 



l 2,R 



it) 



l 2,R 



dt'h R (t-t')\A(t')\ 2 (79) 



where A is the Gaussian pulse and h R is the normal- 
ized Raman response, Eq. (9). The peak value at t = 
will measure the Kerr-like SPM nonlinearity picked up 
through the Raman interaction. Obviously if A is a long 
pulse, it will be 5-like in frequency domain, and the peak 
value will simply be n\ R , i.e. the DC value of the Ra- 
man response. This is what we observe in Fig. 7 for 1000 
fs and even 200 fs input (top two plots): the convolu- 
tion is almost identical to the long-pulse limit where 
n i ~~ * n 2 iil^-W| 2 i which is indicated with a black 



dashed line. When 100 fs or 50 fs pulses are used, the 
bandwidth of the excitation is large enough to couple 
into Raman modes, and the ringing in the convolution 
starts to appear. Note three things: (a) some Kerr-like 
nonlinearity is "picked" up from the Raman term even 
with short pulses, and this must be taken into account if 
short pulses are used for Z-scan measurements, (b) The 
explanation used by Wang et al. that they only observe 
instantaneous nonlinearitics is not quite true. What they 
measure is the electronic plus DC Raman component, 
the former being instantaneous and the latter actually 
some kind of steady state, (c) The fact that they did not 
observe any delayed response on the nonlinearity was ex- 
actly because they did not use a short enough pulse to 
excite the Raman modes. As before we can correct for 
the cascading contributions, which at 520 nm turns out 
to be n^.casc = — 13 x 10~ 20 m 2 /W. Thus we conclude 
that 



n 2,tot 



52 x IQ- 20 m7W 



(80) 



cV, or 60% of the bandgap value in LN. In this regime 
a strong 2PA enhancement of the electronic nonlinearity 
is predicted by the two-band model (the predicted peak 
lies just beyond half the bandgap energy [40, Fig. 1]). The 
relative error on this value is 15%, mainly determined by 
the error of the Z-scan measurement. 

(3) 

Presumably the only measurement of the Xzzzz e ^ cc ~ 
tronic nonlinearity was performed by Wang et al. [41], 
where they used a 0.06% Fe doped MgO:LN crystal 
and a 200 fs 1 kHz Z-scan setup at 520 nm. Fc-dopcd 
crystals are used for enhancing photorefractive effects, 
but on the other hand on a sub-ps scale with kHz rep- 
etition rates they should not kick in. They measured 

Xzzzz = 4 - 96 x 10_13 esu > Le - X (3) = 0.69 x 10 4 pm 2 /V 2 
or «2 = 39 x 10~ 20 m 2 /W. Whether we can assume that 
a 200 fs pulse does not accumulate any contribution to 
the nonlinear refractive index is an open question: they 
claimed that a time-resolved nonlinearity measurement 
did not reveal any non-instantaneous response, and thcrc- 



A. Evaluation of the data 

We choose to use the Raman modes found by Kaminow 
and Johnston for our simulations. Their experiment was 
probably carried out using an sLN crystal, see more be- 
low. One aspect that is relevant here is that we usually 
use cLN and often cLN with 5%MgO doping. The cLN 
crystals turn out to have broader lines than sLN [35, 42], 
in particular when using the 5%MgO doped cLN. Also 
the resonance frequency is shifted slightly [42], which is 
of little importance for our purpose. We therefore use 
the following line widths and resonance frequencies for a 
5%MgO:cLN at room temperature 

v a = {251, 273, 331, 631} cm" 1 (81) 

Y a = {32,16,16,30} cm" 1 (82) 

ri l0 . = {21.2, 19.5, 16.0, 8.4} fs (83) 

T- 2 ,a = {332, 664, 664, 354} fs (84) 
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FIG. 8. Raman model from Kaminow and Johnston [5], i.e. 
Fig. 1 normalized to unity at DC. We modified the line 
widths for 5% MgO:cLN (top) and cLN (middle). The bot- 
tom plot compares the absolute values of the original model, 
presumably sLN (thin dashed) and our modified model for 
5% MgOxLN (thick) and cLN (thin). 



The third mode line width is a bit unsure as we did not 
find accurate information about it for this particular crys- 
tal, so we assumed it to be twice as big as in sLN. This is 
of little importance since it is not such a relevant mode. 



In the same way for cLN we get 

v a = {251,275,331,631} cm" 1 
T a = {28,14,16,27} cuT 1 
r hc = {21.2,19.3,16.0,8.4} fs 
T 2 ,a = {379, 758, 664, 393} fs 



These choices do not change (substantially) the fractional 
Raman values, which are given by Eq. (40), and we repeat 
them for convenience 



(85) 
(86) 
(87) 



f Ra = {0.635,0.105,0.020,0.240} 



(89) 



What does change, though, is the peak-to-DC ratio 
(clearly seen in Fig. 8 bottom): the DC level is found 
using the analysis above, but when we increase the line 
width the peak of the Raman response decreases since 
the area is conserved. At the same time the DC level 
is unchanged, as discussed before. Thus, the ratio will 
increase with an increased line width, and we therefore 
find 



(3) 
Xr' 



max [Xi?,l(peak)] 



(3) 

Xr 



max [Xi?,l(peak)] 

(3) 

Xr 

max [xr^ (peak)] 



= 0.119, sLN (90) 
= 0.17, cLN (91) 
= 0.19, 5%MgO:cLN (92) 



where Eq. (90) is the value from Eq. (43), repeated here 
for context. The value in Eq. (91) is quite similar to the 
values of Schick et al. (where a cLN crystal was used) 
and Phillips et al. (who presumably use cLN crystals). 
We plot the Raman response function in Fig. 8, and the 
bottom plot compares the original model (dashed) with 
the modified models (full). 

We see the line widths are crucial for determining the 
key parameters of the model, and changing them leads 
to different parameters than those used to determine the 
coefficients: the gain measurements determine the peak 
strengths, which in turn through the line widths and res- 
onance frequencies give the DC values for each mode. 
Whether the cubic nonlinear strength changes with crys- 
tal composition is not known, although indications from 
the measurement of the quadratic nonlinearities indicate 
that variations are minimal [38]. However, in order to 

(3) 

extract the DC parameter Xr from a single measure- 
ment of a Raman peak gain, we need to relate its value 
to the DC value using the calculated ratio for that par- 
ticular Raman model. To do this consistently it is impor- 
tant to know which kind of LN the gain measurement is 
performed on. Alas, this information is consistently left 
out. However, educated guesses based on line widths of 
the lowest mode at 251 cm -1 [35, Fig. 3] are that (a) 
Kaminow and Johnston as well as the measurements in 
Boyd's book use an sLN crystal, because the mode has 
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a line width of 23 cm 1 (b) Chunaev et al. use a cLN, 
because the mode has a line width of 28 cm -1 . 

There are measurements of Raman gain and Z-scan 
measurements performed at the same (or similar) wave- 
lengths, namely Eq. (52) and Eq. (48) for the Raman 
gain, and Eq. (77) for the if-scan measurement, which 
gave n{ tot = 38.6 ± 3.2 x 1CT 20 m 2 /W at A = 1.064 /im. 
We now use this value as a reference point for the total 
nonlinearity at this wavelength. The problem lies now 
in determining the DC Raman value from the peak gain 
values. If we use the ratio 0.17 for cLN, cf. Eq. (91), then 
we get for the Chunaev et al. data at A = 1.047 /im [17] 

X ( r = 0.61 ± 0.092 x 10 4 pm 2 /V 2 (93) 
n 2.R = 37.0 ± 6.4 x 10~ 20 m 2 /W (94) 
f R = 0.96 ±0.18 (95) 

Frequency scaling between A = 1.064 /im and A = 
1.047 /xm is considered insignificant. In the uncertainty 
calculations we estimated the error on the Raman gain 
measurement to be ±15%, which is quite conservative. 
Such a high Raman fraction, close to unity, seems quite 
unphysical. However, strong nuclear contributions of 
around 60% to the Kerr nonlinearity has been observed in 
niobium oxide glasses [43], so the result more than any- 
thing is an indication that LN has very strong Raman 
nonlincarities. We also remark the large uncertainty, and 
within two standard deviations one would have a value 
of around f R = 0.6, i.e. along the lines what Schiek et al. 
obtained, cf. Eq. (5). If we use a smaller ratio than 0.17 
then the f R value will be smaller. 

The measurement by DeSalvo et al. leading to Eq. (78), 
n itot = 93 x 10~ 20 m 2 /W, lies close in pump wavelength 
to the Raman measurements of Johnston and Kaminow, 
and we can therefore make a cautious calculation of the 
ratio between the total, Eq. (78), and the Raman DC 
nonlinearity, Eq. (39), to get 

f R = 0.37 ± 0.08 (96) 

This value is therefore a quite good estimate at 532 nm, 
under the simple assumption of a negligible change in 
the Raman nonlinearity (based on the scaling discussed 
above, it should be less than 10%, and the fraction should 



therefore be accurate within this range). This value is 
considerably smaller than the value obtained above, but 
remember that at 532 nm the 2BM predicts a strong 
enhancement of the electronic nonlinearity due to onset 
of two-photon absorption. Since the Raman part does not 
follow this scaling, the result is a much reduced Raman 
fraction. 

In Rcf. [3] we used a 5%MgO:cLN and found that the 
following parameters gave a good agreement with the ex- 
perimental results at Ai = 1.3 /im 

< tot - 45 x 10 - 20 m 2 /W, f R = 0.50 (97) 

There we used the 4-line Raman model that we pre- 
sented in Eqs. (81)-(84). In another recent publication 
[44] we used the same Raman model and parameters but 
a slightly lower ri2 tot value. The value f R = 50% was 
in part inspired by the analysis presented above, where 
evidently very large Raman fractions are possible, and 
in part by looking for agreement between the numerical 
model and selected experimental data. 

Phillips ct al. [2] found a different value at A = 
1.043 /im. They also used n^ tot = 38.6 x 10~ 20 m 2 /W, 
which was estimated in the same way as we did for Eq. 
(77), and now using the previously calculated value Eq. 
(56) we get 

f R = 0.14 (98) 

i.e. a much lower total Raman fraction; this is due to the 
very low peak Raman susceptibility they chose based on 
data from simulations and comparison with experiments. 

A final way of scaling the data is to take a single Z- 
scan measurement, and pair it with all the possible Ra- 
man data. This will give various results for the electronic 
value. Now that we know the electronic value we can 
argue that it scales as cither (a) Miller's delta or (b) as 
the shape of the 2-band model (i.e. we assume that the 2- 
band model can predict the shape if not the exact value) . 
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